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Abstract 

Helical edge states exist in the mixed spin-singlet and spin-triplet phase of a noncentrosymmetric 
(NCS) superconductor. In this article we have considered a planar ferromagnetic metal/NCS 
superconductor tunnel junction and have studied the effect of these helical edge states which 
manifests itself through the charge and spin tunneling conductance across the junction. We have 
shown the behavior of conductances for the entire range of variation of 7 = A_/A + where A± are 
the order parameters in the positive and negative helicity bands of the NCS superconductor. There 
exists a competition between the Rashba parameter a and the exchange energy E ex which is crucial 
for determining the variation of the conductance with the applied bias voltage across the junction. 
We have found a nonzero spin current across the junction which appears due to the exchange 
energy in the Ferromagnet and modulates with the bias voltage. It also changes its profile when 
the strength of the exchange energy is varied. 

PACS numbers: 74.70.-b, 74.45. +c, 74.20.Rp, 75.70.Cn, 75.76.+j 
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I. INTRODUCTION 



Study of the tunneling conductances between different materials coupled together is an 
interesting area of research in condensed matter physics which is also helpful in predicting 
many properties of the systems which are coupled. Ferro magnet ism and superconductivity 
show antagonistic behavior yet when they are coupled together they show some interesting 
transport properties. In the ferromagnetic system the presence of exchange energy creates a 
population imbalance between the two spin sub-bands. When a ferromagnetic metal is cou- 
pled with a normal s-wave superconductor then there occur some interesting phenomenon 
which were absent if the ferromagnet was replaced by the normal metal. In normal metal 
s-wave superconductor junction for a bias voltage less than the superconducting order pa- 
rameter there occurs a process which is called Andreev reflection. By this process the 
incident electron with spin o is reflected back as a hole with same energy but with opposite 
spin a. The change in momentum due to this process is negligible as the Fermi energy 
in normal metal is huge. In case of ferromagnetic metal the presence of exchange energy 
induces an extra momentum change of the reflected hole which causes a reduction in the 
Andreev reflection amplitude. In other words the presence of spin sub-bands reduces the 
availability of the conducting channel and this reduces the Andreev reflection 1 ' 2 . Later this 
idea has been extended to the case of unconventional superconductor/Ferromagnet junc- 
tion and similar observation of the suppression of the Andreev reflection was found through 
the appearance of dip in the zero biased differential conductance 3 . This fact has also been 
utilized to determine not only the polarization of the ferromagnetic material but also the 
nature of the order parameters of various kind of superconductors. It has been studied to 
determine the order parameter symmetry of the singlet d-wave superconductor 4-9 , triplet p- 
wave chiral superconductor 10-17 , mixed singlet-triplet order parameter superconductor 18 and 
recently also in superconducting Graphene system 19,20 . One peculiar feature of this Andreev 
reflection is that it can produce singlet as well as triplet Cooper pairs in the superconducting 
regions provided the superconductor in the opposite side can sustain both the singlet as well 
as triplet order parameters. As we will explain shortly the NCS superconductors appear 
as an ideal candidate here since the order parameter in them is a combination of singlet as 
well as triplet order parameters. This kind of study has been performed by Samokhin et. al 18 
already by calculating the bulk conductance across a ferromagnetic metal/NCS supercon- 



2 



ductor junction. 

The discovery 21 of edge states with definite topological properties in the noncentrosym- 
metric superconductors have intensified the research in this field. The main feature of this 
class of materials is that in them the inversion symmetry is broken. This broken inversion 
symmetry gives rise to some special kind of spin-orbit coupling which is known as Rashba 
spin-orbit coupling(RSOC).The presence of RSOC breaks the degenerate electron band and 
gives rise to two bands with opposite helicities. As a result of this the superconducting 
order parameter instead of having any definite symmetry is a combination of singlet and 
triplet components. This particular combination gives rise to many strange behaviors in 
many of the physical properties of the system. The presence of the Rashba parameter also 
manifests itself by forming the topological edge states. The topological edge states exists 
provided the amplitude of the triplet order parameter is greater than that of the singlet order 
parameter. 22 . The topological edge states are itself interesting and has profound effect on 
the tunneling conductance 22 ' 23 . Many work has been done on them and there exists a spin 
current across the junction of a normal metal/NCS superconductor due to these topological 
edge states, which appears by the application of a magnetic field external to the sample. 

Another more interesting study will be to study the effect of these edge states, present 
in NCS superconductor, across the ferromagnet/NCS superconductor hetero-structure. 
There already exist works on the tunneling properties of these edge states in NCS 
superconductor. 22 ~ 25 . In this article we study the effect of the helical edge states in the charge 
and spin tunneling conductances across a ferromagnetic metal and NCS superconductor 
junction. The NCS superconductor is in its mixed singlet-triplet phase. We have shown 
the variation of the charge and spin conductances with the applied bias and with respect to 
various parameters which are Rashba parameter a, exchange energy E ex , the ratio between 
the strength of the order parameters of the negative and positive helicity bands which we 
call 7 and the barrier height Z. We have found a number of interesting features in our study. 
As there exists a competition between the Rashba parameter and the exchange energy so 
the conductances change accordingly. Another interesting feature is the existence of a tun- 
neling spin current across the junction which appears due to the presence of the exchange 
energy which breaks g c (4>) ^ g c (~ <t>) symmetry (for angle resolved charge conductance) and 
9s(<j>) 7^ —9s(— ( t ) ) (f° r angle resolved spin conductance) by creating an probability imbal- 
ance between the number of up-spin and down-spin incident electrons. This spin current 
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although smaller than the magnitude of the corresponding charge current yet its value can 
be increased considerably by coupling a ferromagnetic material with higher value of the ex- 
change energy. Not only that this spin current or the spin conductance modulates with the 
variation of the bias voltage across the junction of the ferromagnet/NCS superconductor. 
The direction of this spin current is reversed by creating a reverse population density in the 
spin sub-bands. The effect of the parameters a and Z are also studied qualitatively and the 
results are explained. 

This article is organized as follows. In section II we have derived all the formula and 
described the method of our working. In section III we have found numerically the tunneling 
charge and spin conductances from a Ferromagnetic metal to a NCS superconductor due to 
the edge states in the NCS superconductor. We have shown here our results in the form of 
various diagrams and also explained them. We summarize our results in section IV. 



II. TUNNELING CHARGE AND SPIN CONDUCTANCES 

We begin with the Hamiltonian for an NCS superconductor in which Cooper pairs form 
between the electrons within the same spin-split band: 

H= J2 [ekAickA+(A kA 4 A c t _ kA +h.c.)] , (1) 

k,A=± 

where £kA = 6c + Aa|k| for Rashba SOI 26 , £k — h 2 k 2 /(2m) — /i. Here /i, m, A, k, a, and A kA 
denote chemical potential, mass of an electron, spin-split band index (±), momentum of an 
electron, coupling constant of Rashba SOI given by V so = arj^ ■ a with r/k = yk x — xk y and 
the Pauli matrices cr, and pair potential in band A respectively. We choose k y + ik x -wave 
pair in both the bands, i.e., AkA = A A Ak with Ak = — iexp[— ifa]. This corresponds to 
triplet component of pair potential A T = (d k • cr)ia y with eZ k = 2^j(A+ + A_)r7 k , i.e., the 
amplitude of the triplet component A t = |(A + + A_) and the singlet component of the 
pair potential is A s = iA s a y with amplitude A s = |(A + — A_) (Ref.27). Therefore the 
superconductor is purely triplet with k y + ik x -wave symmetry when A + = A_, purely singlet 
with s-wave symmetry when A_ = — A + , and triplet and singlet components with equal 
amplitude when A_ = 0. Therefore the Hamiltonian (1) in the matrix form 28 read as 



where = £k + V so and A k = A T + A s . The solution of the Hamiltonian (2) in the bulk 
is with the energy eigenvalues values ±a/£+ + A 2 and ±a/£ 2 + A 2 ,, in with the Cooper 
pairing between electrons within the same spin-split band. Correspondingly, there are two 
Fermi surfaces with Fermi momenta kp = ^fma/h 2 + a/ (ma/ h 2 ) 2 + 2mfi/h 2 , i.e., kp < kp. 

Now we consider a planar junction between a ferromagnetic metal(a; < 0) and NCS 
superconductor (a; > 0) in such a way that the barrier is along Y direction(x=0). The 
interface is characterized by a potential U(r) = U 5(x) with U as the strength of the 
barrier. The Hamiltonian describing the ferromagnet is given as 

H FM = (-7^ + U(r) - E f m)a - ah (3) 

where a — ± for different spin orientations. This second term in the Ferromagnetic Hamilto- 
nian is called exchange energy interaction and play a important role in the tunneling process. 
Here Ef m describe the Fermi energy in the Ferromagnetic side. We neglect any Fermi energy 
mismatch between two sides of the junction ie. assume that the Fermi energy on both sides 
of the junction are at the same level. We now try to construct the quasi-particle wave func- 
tions in the Ferromagnetic side. To do this we follow the way of writing the wave function 
as described by 18 Samokhin et al. We make further simplification to them by taking E = 
such that the wave vector for particles and holes with any particular a are identical. We 
describe the quasi-particle wave vectors for a =f, J, as, 

fc| = ^2m[E fm + h ] (4) 
k\ = ^2m[E fm -h ] (5) 

An incoming spin incident from the FM side on the interface may be normally reflected 
or Andreev reflected. If we denote by and the normal reflection coefficients , by 
and the Andreev reflection coefficients for an incoming spin a then we can write the 
spin-up and spin-down quasi-particle wave functions as , 
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where we have kept room for all possible spin flip processes at the interface. Here 9^ 
and 6*J are respectively the angles made by the wave-vectors k^ and fcj respectively with the 
interface normal. Similarly the wave function for spin-down quasi-particle is written as, 
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We can also write the quasi-particle wave functions in the NCS superconductor side which 
describe the edge states. They are 
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(7) 



^s(x, y) = e ikyV [e~ K+x {ci ip+e ik F* x + c 2 ^e~ ik ^ x } 
+e- K - x {d l ^e ik '^ x + d 2 ^e- ik r* x }} , 



(8) 



where Fermi momenta along x-direction in two spin-split bands are kp x = yj kp 2 — k^ . 
Quasiparticle and quasihole wave functions 23 in two spin-split bands (±) are given by 
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where ^ = (E - iT + )/A+, ^ = (E - ir_)/A_, and r± = y/&± - £ 2 for an edge state 
with energy E 1 , and sin(0± CT ) = k^/kp, where a stands for the spins of the incoming electron. 
Here ci, c 2 , c?i, and d 2 are the corresponding weights at which these four quasi-particle and 
quasihole states mix, and k± = mT±/kp x are the inverse of the length scales of localized 
edge state for two spin-split bands. 

We have to determine the above coefficients by matching the boundary conditions suit- 
ably. In doing this we have to consider the spin of the incoming electrons. The reflection 
amplitudes can be found out by matching the wave functions and the velocity flux at x = 0: 
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From the conservation of the parallel component of the wave vector one can write for an 
incoming up-spin particle, 

k^sinO = fcjsinflj = k^sinO^ = k F+ sin6 + ^ = k F _sin6_^ (13) 



Similarly for an incoming down-spin particle we have, 



k e ^sin9 = k^sinO^ = k^sinO^ = kp+sin9 + i = kp-sinO-i (14) 

This is a peculiarity that appears in the case of a Ferromagnetic metal is that the presence 
of exchange energy causes an imbalance between the density of states of the up and down 
spin electrons. As a result the probability that an incoming electron is spin up(Pj) differs 
from that of an electron with spin down (Pi). This probability factor is given as P a = 
l(l + ah /E FM ). 

Extending the generalized BTK formalism 29 in our case we can write the angle resolved 
charge (G c f(-E, 9)) and the spin(G s i(E,0)) conductance as(for an j electron incoming) 

G c i(E, 6) = 1 + (IrJl 2 - W?) + (^)d^l 2 - ia 2 ) (15) 

G.i(E,6) = (|rf| 2 - + (^)(l^| 2 - |r^| 2 ) (16) 

Similarly the angle resolved charge(G c |(-E, 9)) and the sp'm(G s i(E, 9)) conductance as(for 
an I electron incoming) 

G cl (E,9) = 1 + (|^| 2 - |r e U| 2 ) + (^)(|ri T | 2 - |rf| 2 ) (17) 

G si (E,9) = (|rf| 2 - |r e ^| 2 ) + (^)(l^| 2 - |rf | 2 ) (18) 

Here the factor for f-spin case and for j-spin case is the fraction by which the 
Andreev reflection process and the spin flip process at the interface get suppressed. We can 
find out the angle integrated values of the above quantities by integrating over the angle of 
incidence. We defined them as, 

G aT {E) = -±-l d9cos9G ca (E,9) (19) 
Gn Je c 

G sa (E) = -J- / d9cos9G sa (E,9) (20) 
Je c 

where Gn is the tunneling charge conductance from normal metal to normal metal and 9 C 
is the critical angle of incidence. We define two critical angles of incidence here. One is 
9 C which comes from Equ (13) as 9 C = sin~ l (-±). For an incident angle 9 > 9 C there is no 



Andreev reflection process . There is another critical angle which we denoted by 9 CS and 
defined as 9 CS = sin~ l (^-). For an angle of incidence greater than this angle there will be 
no transmitted particle. This restrictions on 9 comes from the situation when the incoming 
electron is with up-spin. It is also noticeable that similar restriction does not arise for a 
down-spin incoming particle. In fact we can integrate the tunneling conductance for down 
spin particle over the entire angle of incidence, but since the up-spin particles impose some 
restrictions on the angle of incidence we have to follow that restriction in case of down-spin 
also. So the actual critical angle will be the smaller one between 9 C and 9 CS . To determine 
this the relative magnitude of a and E ex is important. It can be easily shown that for 
a < E ex , 9 C < 9 CS . So the actual critical angle should then be 9 C . In the reverse case 9 CS 
would define the appropriate limit of integration. 

Finally we define the total charge and spin conductance as, 

G C {E) = ^ a P a G ca {E) (21) 
G S {E) = ^ a aP a G sa (E) (22) 

where a is positive for up-spin and negative for down-spin incoming electrons. 

III. RESULT AND DISCUSSION 

(A)Role of Andreev reflection in the tunneling process 

In the Andreev reflection process the energy as well as spin is conserved. When a normal 
metal is coupled to a s-wave superconductor then an electron with up-spin after Andreev 
reflection returns as a down-spin hole. As a result a Cooper pair is formed with charge 2e 
and spin zero ie. in singlet state. The presence of exchange energy in ferromagnetic metals 
gives rise to two different spin sub-bands with different population. The noncentrosymmetric 
superconductor can sustain both singlet as well triplet order parameters due to the presence 
of Rashba spin-orbit coupling and two bands are of different helicities. As a result when 
these two materials are coupled together then on the ferromagnetic side the incident electron 
with spin a coming from the a sub-band may be Andreev reflected as a hole in the a = —a 
sub-band with spin a = —a or as a hole in the a sub-band with spin o. In the former 
process there is a tendency of formation of singlet Cooper pairs with opposite spin but in 
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FIG. 1: (Color online) The variation of charge conductance g c and spin conductance g s (normalized 
with Gjsf) with incident angle 9 for different values of exchange energy E ex .The parameters Z = 
5,7 = 0A,a/vFs = 0.1, eV = 0.2A+ are chosen. The exchange energy are 0.0 (black solid line), 0.2 
(red short-dashed line), 0.4 (green dashed-dot line). All in the unit of Ef m . In (A) the population 
in the up-spin sub-band is greater than the down-spin sub-band and in (B) the population in the 
down-spin sub-band is greater than the up-spin sub-band. 

the later process triplet Cooper pairs with spin (a, a) is formed to conserve the spin in the 
Andreev process. Although instead of propagating as separate entity they ultimately form 
quasiparticles which is a combination of both of them. So in this process some spin current 
is transported due to the edge states present in the noncentrosymmetric superconductor. 
Since in our sample the spin-up band is given some preferential population over the spin- 
down band so the triplet pair with both spin down will be produced preferentially. As a 
result a negative spin current will flow as is seen from the asymmetry of the Fig. 1 which 
is opposite to the charge current. This spin current can be reversed by creating a greater 
population in the down-spin band. Although in other figures we have plotted the absolute 
value of this negative spin current so they appears positive. 
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(B)Discussions of the numerical results 

We have numerically solved equations 11,12 to get various coefficients and later they 
were numerically integrated to find out the values of the various angle integrated quantities. 
In Fig. 1 we have shown the variation of both spin and charge conductances (normalized 
with Gn) with the variation of the angle of incidence 9. Since they are the values of the 
angle resolved conductances so we denote them by g c and g s . The presence of the exchange 
interaction is incorporated in such a way that it creates a preferential population to the spin- 
up sub-band. Also there is an preferential probability for the incident up-spins. This extra 
probability factor in ferromagnetic metal makes an 6 dependent asymmetry in the charge 
as well as in the spin conductances and breaks the symmetry g c {(f>) ^ g c (—<f>) (f° r angle 
resolved charge conductance) and g s (<f>) ^ —9s(—<fi) (for angle resolved spin conductance) 
. As a result a spin current appears along with the conventional charge current across 
the junction due to the presence of the edge states in the NCS superconductor. This spin 
conductance arising due to the presence of the edge states and without the application of 
the magnetic field also modulates with the variation of the bias voltage across the junction. 
With the increase of the exchange energy the probability for up-spin incidence increases and 
the spin conductance looses the symmetry more and more as is seen in Fig. 1. We have 
shown in these diagrams the variation for the case when exchange energy is absent. In this 
case the variation is odd with respect to 9 and the net spin current is zero when integrated 
over the angle of incidence. We have also shown the result (figure (B)) when the population 
in the down-spin sub-band is greater than the up-spin sub-band. 

In Fig. 2 we have plotted the variation of angle integrated charge and spin conductances 
along with the bias voltage, normalized by the normal metal/normal metal tunneling con- 
ductance G^. The structure has almost the same dome shaped feature as in the case of 
normal metal/NCS superconductor 23 . Here in each plot we have shown the variation of the 
charge and spin conductances for different values of the exchange energy. As expected the 
charge conductance decreases along with the increase of the exchange energy E ex . The rea- 
son for this is that with the increase of the exchange energy the energy difference between the 
two opposite spin bands increases and as a result the probability of both normal as well as 
Andreev reflection decreases and hence there is a decrease in the charge conductance. Also 
with the increase of the exchange energy the critical angle of incidence 6 C decreases along 
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FIG. 2: (Color online) The variation of charge conductance (left panel) and spin conductance (right 
panel) with bias energy eV for different values of exchange energy E ex . We have also shown the 
variation for different 7. The parameters Z = 3 and ajvp = 0.1 are chosen. The exchange energy 
are 0.2 (black solid line), 0.4 (red long-dashed line), 0.6 (green short-dashed line). All in the unit 

with the critical angle inside the NCS superconductor 8 CS .However the spin conductance has 
somewhat opposite trends. Initially with the increase of the exchange energy as the energy 
difference between two opposite spin sub-bands increases the probability that an electron 
with a certain spin returns with the opposite hole after the Andreev reflection is reduced. 
Also with the increase of the exchange energy the probability that the incoming particle 
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FIG. 3: (Color online) The variation of charge conductance G c and spin conductance G s with bias 
energy eV for different values of Rashba parameter a/vF s . The parameters Z = 3,7 = 0.3 and 
E ex = 0.2Ef m are chosen. 

is with spin up increases. This seriously affects the spin conductance but have little effect 
on charge conductance. As a result the formation of Cooper pairs with the parallel spin 
triplet pairing becomes more viable compared with the singlet pairing and this gives rise to 
an increased spin current along with the increase of the exchange energy. This spin current 
though rapidly falls down when the exchange energy is increased more after a certain value. 

We have also shown the feature of the charge as well as spin conductance for different 
ratios of the 7 = For each 7 a dip appears in the conductance diagram exactly at the 
position where the bias voltage eV equals 7. In fact this feature is universal and always 
appears whenever we will plot charge and spin conductances with the variation of eV for any 
value of 7. One interesting feature of the figures for 7 7^ 1 is that the variation of the charge 
conductance for eV > 7 changes more with the increases of the exchange energy compared 
to that within the dome shaped region. It happens due to the fact that quasi-particles with 
eV > 7 in the NCS superconducting region effectively see one band, the bound states of the 
other band simply becomes inaccessible. Within the dome shaped region however both the 
bands take active part in the conduction process and as a result the variation of the heights 
of the peaks of the charge conductance for small variation of the exchange energy becomes 
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FIG. 4: (Color online) The variation of charge conductance with bias energy eV for different values 
of barrier strength Z. 

negligibly small. 

In Fig. 3 we have shown the variation of the charge and spin conductances for different 
values of the Rashba parameter a keeping the other parameters fixed. This figure is also 
interesting. This figure shows that with the increase of a both charge and the spin conduc- 
tance decrease. The reason is obvious. With the increase of a the critical angle 9 CS decreases. 
We would like to mention here that for a > E ex the angle 8 C exceeds 8 CS . As a result we 
have to perform the integration between±# cs . So there is a competition between the Rashba 
parameter a and the exchange energy E ex . When a < E ex then 9 C is the critical angle but 
whenever we have the reverse situation i.e. a > E ex then 8 CS becomes the critical angle 
for integration. Finally in Fig. 4 we have shown the variation of the charge conductance 
with the variation of the barrier height. The conductance decreases with the increase of the 
barrier height and the behavior is as usual. 
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IV. SUMMARY 



We summarize by saying that the helical edge states 22 exist in a noncentrosymmetric 
superconductor provided the triplet-pair amplitude is larger than the singlet-pair amplitude, 
i.e., when < 7 < 1. We have studied the effect of these edge states on the tunneling charge 
and spin conductance between a ferromagnetic metal and a NCS superconductor. The 
presence of the exchange energy breaks g c {4>) 7^ g c {—<t>) symmetry (for angle resolved charge 
conductance) and g s (4>) 7^ —9s(—<fi) (for angle resolved spin conductance) by creating an 
probability imbalance between the number of up-spin and down-spin incident electrons and 
as a result there exist a tunneling spin current along with the usual charge current due to the 
presence of the edge states in NCS superconductor. Not only that, this spin current can be 
tuned by tuning the bias voltage across the junction and it changes shape with the change 
in the exchange energy in the ferromagnetic metal. The direction of this spin current is 
reversed by creating a reverse population density in the spin sub-bands. We have shown the 
variation of both the conductances with the different values of the ratios between the two 
order parameters which we call 7. The variation of the charge as well spin conductances with 
the change in the Rashba parameter a and the barrier height Z are also presented. From 
these plots we can conclude that there is a competition between the Rashba parameter a and 
the exchange energy E ex which manifests itself through the presence of two critical angles 
6 C and 9 CS and the variation of charge and spin conductances . The plot of the variation of 
conductance with Z shows no unusual behavior. 
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